The time evolution of permutations under 

random stirring 



Balint Veto 

Abstract. We consider permutations of {1, ... , n} obtained by [y^t\ indepen- 
dent applications of random stirring. In each step the same marked stirring 
element is transposed with probability 1/n with any one of the n elements. 
Normalizing by we describe the asymptotic distribution of the cycle struc- 
ture of these permutations, for all t > 0, as n ^ oo. 

1. Introduction 

We consider the following random stirring mechanism: n numbered balls are given in 
the beginning on their corresponding numbered places. In each step, independently, 
the first ball, which is referred to as the stirring particle or stirring element, changes 
place with one of the n balls or stays unchanged with probability 1/n. We investigate 
that permutation which brings the balls from their initial place to their place after 
i steps. 

Formally, let tt'"^ (i) = T^*'"' o T^^\ o • • • o T^*-"-* be a permutation acting on the 

set [n] :— {!,..., n}. The permutations (T-''"')^^ are chosen independently with 
uniform distribution from the n — 1 transpositions moving the stirring particle and 
the identity permutation. 

Let cr be a permutation of a finite set S, i.e. an 5 — > S* bijective function. The 
cycles (orbits) of a are the sets of form {v, <j{v), o'^(w), . . . } C S for some v (1 S. 
The set S is the disjoint union of its cycles. The cycle structure of cr is the sequence 
of the cardinalities of the different cycles in non-increasing order. 

In our case one of the cycles can be distinguished from the others (namely 
the cycle of the stirring element), which will be called the active cycle. For the 
total description it is enough to determine the distribution of the cycle structure 
of the permutation 7r*^")(i) (regarding the active cycle separately). This gives the 
distribution of the conjugacy class of tt^"^ (i) restricting ourself to the conjugation 
with permutations fixing the stirring particle. The distribution of tt^") (j) is uniform 
within a fixed conjugacy class. 

We encode the permutation 7r(")(i) with the vector C^^^i) := {C^''\i) , c["'\i) , 
C2"\i), . . . ) where C'^\i) denotes the length of the active cycle, C^\i), C'^\i), . . . 
the lengths of those cycles in non-increasing order which are already moved by one 
of the transpositions (Tj-"^)*^j^. Other cj"''(i)-s are 0. (C^"' (i)),^g is a process on 
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Figure 1: The metric on S 



the state space 

S := {{so,si, S2, ■■■)■■ Sn eR, s„ > n = 0,1,2,..., 

si > S2 > • • • > s„ > . . . and sj > for finitely many j} 

with the distance 



d{A, B) := sup ■ 



fc = 0,l,2,. 



(1) 



(2) 



where A = (Aq, Ai, A2, . . .) and B — {Bq, Bi, B2, . ■ .) are elements of S. (See Figure 
nj) The ranking is not a natural part of the problem, but it facilitates studying the 
model. 

At each step after applying a random transposition two types of changes may 
happen in the cycle structure: merging of two distinct cycles or splitting of a cycle in 
two. While different transpositions are applied (meaning that the stirring 

particle chooses a new element in each step until i), the cycle decomposition of 
7r*^"^(i) contains only fixed points and the active cycle, which increases by one in 
each step: C'"^(i) = (i + 1, 0, 0, . . . ). If a transposition recurs, then the cycle splits 
in two, one of which will be the new active cycle. If there are already more than 
one non-trivial cycles in the decomposition, then the active cycle can merge another 
cycle. (See Figure [2) The model realizes a coagulation-fragmentation process. 

A reduction of the problem is to study the coagulation and fragmentation events 
of the cycles together, because both of these events happen when the stirring element 
steps to a place already visited. We investigate this simpler question first. Then we 
introduce a continuous time process on S, which turns out to be the limit process. 
The convergence is proved by coupling. In Section 4. we show that the station- 
ary distribution of the underlying split-and-merge transformation is the adequate 
modification of the Poisson-Dirichlet distribution. (See the definition later.) 

A similar model is studied by Schramm in [Q. He chooses {Ti)°li to be inde- 
pendent random transpositions with uniform distribution from all possible transpo- 
sitions of the set [n] . The limit distribution of the proportions of the giant cycles 
in the permutation n(t) = Tt o ■ ■ ■ o Ti after t — cn steps as n ^ 00 is identified 
(where c > 1/2 is a constant). 

This result is in accordance with the classical theory of the random graphs 
derived from Erdos Let us consider the random graph G{t) on the vertex set 
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Figure 2: Coagulation and fragmentation of cycles 

[n] where {u, v} is an edge in it if and only if the transposition (u, v) appears in 
{Ti, . . . , Tt}. By the Erdos-Renyi Theorem jZj the graph G{t) has a giant connected 
component only in the case t/n — c > 1/2 similarly to the condition on the random 
permutations. (For random graphs and random graph processes see jH] and [T2]l. 

In Schramm's paper the vector of the cycle sizes of II{t) in non-increasing order 
normalized by the magnitude of the giant connected component of G(t) converges in 
distribution to Poisson-Dirichlet distribution with parameter 1 after t = cn steps 
(c > 1/2) as n ^ cx). That is the limit distribution of the relative cycle sizes in 
a random permutation chosen uniformly from all permutations of [n] as n —>■ oo. 
Thus for large n the permutation n(t) behaves on the giant connected component 
of the Erdos-Renyi graph G{t) as a uniform permutation. 

Our paper is motivated by Toth in connection with the quantum-physical appli- 
cations of the problem ^^l- Angel analysed Toth's random walk model on regular 
trees in For similar random stirring models see also [2] and 0. 

2. Return times of the stirring particle 

The movement of the stirring particle is a random walk (i?|"'')^o ^he set [n], 
which is homogeneous in space and time. Let 

y/") := #{fc : fc < z, 3j<k: = B^^^} (3) 
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be the number of the returns until the zth step to places already visited by the 
random walk (-B]"^)!^o- ^^^^ include those steps when the stirring particle 
keeps its place. 

After the ith step the stirring element has already visited exactly i + 1 — T^*-"-* 
places (including the starting point) , so the transition probabilities of the Markov- 
chain (V^/"^)^o 

■ _L 1 T/^^^ 

p (vl^l - y/") - i\vt^) = 1 - p (v^l - y/"' = o\vt^) = \^ ' ■ (4) 

In order to get a non-trivial limit distribution the time of the processes should 
be accelerated. As opposed to Schramm pi], in Theorem 1 the scaling is ^/n. This 
means that we describe the beginning of the evolution, because after 0{^/n) steps 
the bulk of the elements is still unchanged. Simultaneously we normalize the cycle 
sizes with and we let n ^ oo. 

From now on we investigate the limit of the vectors f *^' '(1V"*J) A as n ^ oo, 
* V V" /t>o 

where the division is meant coordinatewise, namely — — [y^*^"^ {— — ~7^~^i 

— — • ■ • )• Elementary calculations, similar to the classical birthday problem, 
give the following Hmit distribution of the returns. For limit theorems related to 
generahzations of the birthday problem see also fSj. 

Proposition 1. Let (Vt)t>o be an inhomogeneous Poisson point process with intensity 
p{t) = t. Then 

{vl%,^)t>,A{Vt)t>, {n~.oo) (5) 
in terms of the finite dimensional marginal distributions. 



3. Coupling 

Much more can be stated for the above model. Not only (V^^^^^j )t>o, but the se- 
quence of the processes — ^7=-^^^ converges. Moreover by means of coupling 

a stronger type of convergence is realized. 

The limit process is a natural continuous extension of the discrete processes 
(C("^(i))^Q. For large n the active coordinate c'^\i) increases in the bulk of the 
steps (when no split or merge occurs). In the times of jumps of (V^''"'')^o ^ ^P^i^ 
or a merge happens depending on the proportions of the cycle sizes as follows. The 
probability of a split in the ith step, conditionally given that the stirring particle 
returns to a place already visited, is 

(6) 
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The conditional probability of the merge of the jth cycle and the active one is 



We define an S valued continuous time stochastic process C(t) = {Co{t),Ci{t), 
C2{t),...) with cadlag paths, which imitates the above process. It is built on a 
Poisson point process (14)t>o with intensity p(t) = t. Similarly to the discrete 
processes (C^"''(i))^g at the times of jumps of {Vt)t>o a split or a merge event 
occurs with probability proportional to the coordinates of C. 

The initial state is C(0) := (0,0,0,...). The evolution of the process is the 
following: the coordinate Co{t) increases with constant speed 1 between the jumps 
of {Vt)t>Q- Let Tfc be the fcth time of jump of {Vt)t>o, in other words Vt^ = k and 
Vt^- — lirngj^o ^fc-e — k — 1. Let (C/fc)^i be i.i.d. random variables with uniform 
distribution on [0, 1] independent of (Vt)t>o. One of the next two actions occurs at 
time Tfc. 



1. Split: If 



rr ^ Co(rfc-) 

l^m=0 '-'m[Tk~) 



then let Co(rfe) := Uk Zl^^o C'm(rfe-), and the sequence {Cm{Tk))'^=i will be 
the collection of (Cm(Tfe-))^^i and Co(Tfe-) - Uk I]m=o C'm(Tfc-) rearranged 
in decreasing order. 

2. Merge: Otherwise a unique index j > 1 can be chosen a.s. via 

12in=0 CmjTk — ) ^ jj^ ^ X]lt=0 Cm(Tfc — ) 

Let Co(Tfe) := Co(Tfe-) + Cj (r^-), and Cm(rfe) := Cm{Tk-) if 1 < m < j, and 
Cm{Tk) '■— Cm+i{Tk—) li m > j restoring the decreasing order. 

Observe that J2m=o^mit) = t, but we did not use it to simplify the formulas 
(jHJ and (EJ in the above definition because the analogous discrete assertion is not 
true, compare with © and 10. 

The main result of this paper is that the normalized discrete processes converge 
in probability to {C{t))t>o in the following uniform sense in terms of the distance 
defined by 

Theorem 1. There exists a probability space (il,JF, P), on which the discrete pro- 
cesses (C'^"^(i))^Q n — 1,2, . . . and the continuous time process {C{t))t>o can be 
jointly realized so that ifT>0 is fixed and f(n) is any function tending to infinity 
with n, then 



sup d C(t), ^' < ^ ^ 1 

Vo<t<T V y/n J y/n J 



as n ^ oo. (10) 
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3.1. The convergence of the return process 

Let (iljJF, P) be such a probability space where a Poisson point process {Vt)t>o 
with intensity p(t) — t and the i.i.d. random variables (C/fc)^i and with 
uniform distribution on [0, 1] are given independently of each other. 

We have constructed the process {C{t))t>o from (Vt)f>o and {Uk)^^i earlier. We 
first re-create the processes (V^/"^),^o with the appropriate distributions on the new 
probability space (17, J^, P). The main idea of the construction is that we observe 
the process {Vt)t>o in ^ long time intervals. 

Let Xf"^ 1 - VVi > i = 1,2 . . . n = 1, 2, ... be the indicators 

of the increase of the process {Vt)t>o, which are Bernoulli random variables with 
respective parameters 

p!"^ = l-exp(-^)^l + o(4). (11) 
The required parameter for the increase of is 

(n) 



,t' = --^- (12) 

n n 



We define the values of V^'"-* for fixed n with induction on i. Let Vq ""^ :— n ~ 
1,2,... and 




/ (n) _ , _ 

+1 (pf) < gi")) a (x^ = o) a ( zi") < ^YT^rr ) ■ (i^) 

We define y/") -(//."i^ + 

It is easy to see that the distribution of the new (T^/"^)^o accordance with 
. Later on we say that a correction happens if the products of the indicators in 
l(T3ll do not disappear. We will see that the total probability that a correction ever 
occurs is small if n is large enough. This gives an alternative proof of Proposition 
1. 



Lemma 1. Let T > be fixed and denote = ro,ri,...,TK the random times 

(n) (n) (n) 
'^O ' ^1 1 ■ ■ ■ ' \,(") 



of jumps of the process (Vt)o<t<T and denote = rQ"'',r}"'', . . . ,t';"„^) that of the 



discrete process (^[^y^jj )o<t<T defined above. Then for sufficiently large n with 

probability close to 1 the number of the jumps are equal: k = k^"' . Furthermore, 
there exists a bijection between the jumps of the processes in such a way that 

\rk-Ti''^\<4= k^l,...,K (14) 

holds with large probabili 
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For technical convenience we introduce the following events for fixed e,5 > 0: 

E, -.^ {Vj^;^^^ < K, n^N„N, + l,...}, (15) 

where K,. is a sufficiently large constant and is a threshold satisfying F{E^) > 
1 — e. It makes sense by Proposition 1. Let 

Ms := { min {rfe - Tk-i} > 6} n {Vt - Vt^s = 0} (16) 

fc:Tfc<T 

where is the time of the fcth jump of the process (Vt)o<t<T and tq = 0. It is 
elementary that lime^o P(-E^e) = fima^oPl^'^i) = 1- 

Proof of Lemma 1: By (jl(i|l . on the event Ms the increment of the process (Vt)o<t<T 



on any interval 



i i+1 



does not exceed 1 if n > , hence V i — Vi-i — X, 



(«) 



Since v'^'f' — V"f "i = y/"^ , it is enough to prove that 

P({3i < [V^TJ : Xf") ^ i:/"^} n n M5) ^0 (n ^ 00) (17) 

for afi fixed e, (5 > 0. 

On the event xf""* = 1 can be true for at most many indices i. So the 
probability of the correction in the cases p,-"-* > g^^"-* satisfies 



using the power series of the exponential function and the equations ifTTll and lfT2|l 

in) (n) 

estimating p- and ql . If we add this at most many times, then the sum still 
goes to as n ^ oo. A similar calculation shows that for an i, for which p^"-* < q^"^ 
holds, the probability of the correction is at most 

Summing up for i = 1, . . . , the total probability still tends to 0, as required. 



3.2. Splits and merges 

With the processes (V^/"^)^o have determined when a split or a merge occurs, 
our task is now to define how it should happen. Similarly to the definition of 
the limit process (C(t))f>o we can prescribe the evolution of the discrete processes 
(C(")(i))=^o with the use of the same independent uniform random variables {Uk)'^i 
as follows. Let C^"^(0) := 1, C^"^(0) := m = 1, 2, .... The evolution of the 
process C'-"-' in the steps i = 1,2, . . . is described below: 
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• if y/"^ - V}"1 = 0, then C^"^^) := C^''\i - 1) + 1 and other coordinates 
unchanged, 

• if X^/"-* — — 1 and v}'^^ = k, then the uniform random variable Uk 
determines a unique index j with probabiUty 1 as in Q via 

J2m=0 — 1) ^ jj <- J2m=0 — 1) .^q-) 

Similarly to the definition of the limit process 

1. i = 0: split. If UkYlZ=o^m\i — 1) < 1, then let everything be unchanged: 
C(")(j) := C(")(i - 1), let us call this case fictive split (corresponding to 
the event that the stirring particle keeps its place). Otherwise CQ"'\i) ■— 
[Uk Em=o ^ («-l)J ' let the broken fragment C^"^ {i-l)-[Uk E^=o ^ 

1)J add to the collection of nonactive pieces {Cm\i — to form the new 

ranked sequence (Cm'(j))m=i- 

2. j > 0: merge. Let Co"^(«) := cl^\i — 1) + cj"''(« — 1) and for the re-ranking 
C'^\i) := Cl-:\i - 1) if < m < j, and ci"^(i) := - 1) if to > j. 

It is easy to show that this new definition of (C^"^(i))^g provides the same dis- 
tribution as in the model generated by transpositions, so we prove the convergence 
for these processes. 

Proof of Theorem 1: Let e, 5 > be fixed. Let A„ denote the event that the 
assertion of Lemma 1 holds for (V'l^^^^j )o<t<T- We restrict ourselves to the events 
n Ms n An . Let us define a measure (which is not a probability measure) on the 
sets B £ T: 

^e.sAB) ■■= F(B C^E,C^ Ms n An). (21) 

By Lemma 1 it is enough to show that for fixed e, (5 > the processes C(t) and 
— — ^j^^*^"* are sufficiently close to each other for large n except a set with ^e.s.n- 
measure tending to as rt ^ cxd. The proof consists of the following steps: 

1. We estimate the increase of the distance between C{t) and — — ^j^^*^"* between 
two successive split or merge events. 

2. We introduce those events when the distance under discussion cannot be es- 
timated: the awkward events (defined later) and the fictive splits. We show 
that they have small probability. 



3. On the complementer event, which has probability tending to 1 as n ^ oo, we 
show that a 
very much. 



show that a merge does not increase the distance between C{t) and — — ^j^^*^-* 
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4. We do this also for the spHts. 

5. We summarize the estimates. 
Step 1. Let 

(22) 

denote the distance between the discrete and continuous processes before and after 
the time of the fcth split or merge. (Recall that is the time of the fcth jump of 
(Vt)o<t<T and r^""* is that of {V^^^^ )o<t<T, which are close Pe,5,„-almost surely by 
Lemma 1.) 

While no split or merge occurs, the distance between the processes does not 

fn) 

increase very much. From Lemma 1 the difference between and r^; can be at 
most The discrete processes i^—-^^^^)t>o change only in the times which 



are multiples of -y^. Thus Pe,5.„-almost surely 

dfe < d+_, + ^. (23) 



2 

n 

Step 2. From now on we investigate only the split or merge points of the 

-(") \ 



processes. At the fcth time of jump of (Vt)o<t<T and (V,*-"^ , )o<t<T we choose 



with the help of Uk one of the components of C(rfc— ) and C^")([-ynr^."''j — 1) via 
(jnj and lf2(Hl . Let us call the possibility that these components are of different 
indices an awkward event. If an awkward event or a Active split (meaning that 
Uk J2m=o ^™ ^ (* ~ 1) < 1) occurs, then we cannot estimate d{C, C^^^ / ^/n). We will 
see that these events have probability tending to as n ^ oo. 

We can choose the components of C and C(") as follows. We set the coordinates 
of the vector C(r/c— )/ ^^^q Cm(Tfc — ) to the real line from the origin one after 
another, which gives a partition of the unit interval [0, 1]. We do this also with the 
coordinates of 

C(")(L^AIr(")J-l) ^^^^ 



Let Wk denote the set of those points in [0, 1] which are covered by the coordinates 
of C and C^") of different indices. The probability of the awkward events (which is 
an upper estimate for their Pg^^^n-measure) is exactly the Lebesgue measure of Wk- 
We know that X]m=o ^m{t) = t for all t > 0. From the construction 

iv^^\-^^^^r..oCi^\ivEti)^^ if .e[o,n (25) 
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because at the spUt or merge points (occurring at most K,. many times) the total 
length of the discrete process does not increase. From this 



2^ c™(Tfc-) -= 



< 



Tk 



< 



/n X n 



(26) 



follows using Lemma 1. 

From the above it is an elementary exercise to show that the distance between 
the corresponding dividing points of the partitions of [0, 1] generated by C(Tfc— )/Tfc 
and by the vector Ij24|l can be, respectively, at most 

J- , j<-.+3 

Tk 

Since the number of coordinates is at most iCe , this provides the following upper 
bound: 



Leb(VFfc) < 



Tk 

-vOO 



(27) 



where we used the fact that Tk = X]m=o Cra{Tk — ) > S holds for A: = 1, 2, . . . on the 
event Mg. This yields 



IPe.a.n (awkward event at Tk) < Leb(Wfe) < 
Furthermore 



Pe, 5, n (Active split at t^""*) < 



1 



< 



if n is large enough by Ij25|l . So we conclude that 

Pe, (5, n (awkward event or fictive split at the fcth split or merge point) 

K _ Kl + 3K^ 

- ~r"fe ^ ■ 



(28) 



(29) 



(30) 



Step 3. In the case when the random variable Uk chooses the same components 
of C and C^") and it is not the active coordinate, i.e. there is a merge in both 
processes (see Figure El, then 



4 < 3dfe . 



(31) 
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Figure 3: The piece C2 merges Co parallel with the C^"^ - C^"^ coagulation 



X' , 
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B 



B' 
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Figure 4: Split: the broken pieces from the coordinate are X and X' which have 
to be moved to places A' and B 



Step 4. If a (non-fictive) split occurs in the discrete and continuous processes, 
then using inequality ll26|l we have 



CQ{Tk) 



< 



m— 



< 



K, + 3 1 



(32) 



This is why the broken pieces from Co{Tk~) and CQ^\ly^T^f."'' \ — (denoted 
by X and X' on Figure can differ at most + ■ the difference can be -^J^ 

between the left end points and at most c?^ between the right end points. 

It is possible that the two broken pieces do not come to the same place in the 
decreasing order of the coordinates. This case is shown on Figure^] Then we move 
first both X and X' to the closer of the final places of them in the decreasing order 
(to the places A and A' on the figure). Because |yl — A'l < d^, the result is two 
vectors (the modifications of C(Tfc— ) and C^^^^ l^/nrj.!^^ \ — l)/\/n, but one of them is 
not necessarily in decreasing order), which have d{-, •)-distance at most 2d^ + 
more then before this modification. 

In the second step we move X from AtoB (see Figure ^ . The lengths of the 
parts between A and B are at least \X\ and at most \X'\ + 2d'^ < \X\+3d'i^ + 

So any two of these parts have lengths differing at most 3d^ + ■ Swapping X 
always with its right neighbour until hitting place B, the number of the swaps is at 
most Kg, and at each swap the distance can increase at most Sd^ + ^^'^ , so we 
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have 



Step 5. Summing up the estimates ll23| . ll3T| and we get easily the following 
recursive bound: 

d+ < max(3, 2 + 3K,)d+_, + ^= '-^ ^ ad+.i + ^. 

(34) 



Hence supg<j,<^^ d^, < [J2i=oba^j "y^- Considering the results of steps 1 and 2 
the assertion of the theorem follows. 



4. Stationary distribution and generalizations 

It is a natural question to identify the stationary distribution of our stirring process. 
This means that we look at the asymptotic behaviour of the process (C^"-* ([nij )/n)t>o. 
Observe that the time scale is of order n, i.e. the time scale when the stirring element 
has already visited the bulk of the n places. This setup is the same as that of the 
problem studied by Schramm in [n], but different from the phenomenon described 
by Theorem 1. 

In this section we consider the following split-and-merge transformation corre- 
sponding to the stirring generated by random transpositions. Let C — (Co,Ci, 
C2, . . . ) G S be a random probability distribution, i.e. Cm = 1 almost surely. 
Co is the active component. Let U he a random variable with uniform distribution 
of [0, 1] which is independent of C. If J7 < Cq, then the Co splits, i.e. the new active 
component will be U and (Cg — C/, Ci , C2 , . . . ) will be the remaining components af- 
ter restoring the decreasing order. If X)^=o ^™ < U < X)m=o ^^^^ C'o merges 
with Cj similarly to because Cm = 1 • 

In limit theorems of random partitions and permutations the following distri- 
bution appears often. Let the random variables Wi,W2, ... be independent with 
uniform distribution on [0, 1]. Let (Qi, Q2, ■ • ■ ) be the decreasing rearrangement of 
the random variables 

(Pi,P2,...) iWi,{l~Wi)W2,{l-Wi){l-W2)W3,...). 

Then the random sequence (Pi,P2,---) has GEM(l) distribution after Griffiths, 
Engen and McCloskey. {Qi,Q2, ■ ■ ■) has Poisson - Dirichlet distribution with pa- 
rameter 1, abbreviated PD(1). For more about this family of distributions see |^. 

Let (pi,p2,---) be a random probability distribution. We construct its size 
biased permutation. Let J7i, C/2, . . . be i.i.d. uniform random variables on [0, 1] in- 
dependently of (pi,p2, ■ • • )• Let Ij be the unique index for which X)i=i^-Pi — <^ 
Si=i Pi- Let Jk denote the kth smallest integer m satisfying ^ {/i, /2, . . . , Im-i}- 
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Then the vector (p jj , pj^ , • ■ • ) is called the size biased permutation of (pi , p2 , • ■ • ) • 
It is well known that the size biased permutation of a random partition with PD(1) 
distribution has GEM(l) distribution. See also pjQ'l. 

Consider the following probability distribution on S. Let [Qi,Q2, ■ ■ ■) have 
PD(1) distribution. Let Co be a size biased part from [Qi,Q2, ■ ■ ■) (i-e. the first 
component of the size biased permutation of {Qi,Q2, ■ ■ ■)) corresponding to the 
active cycle and the rest (Ci, C2, . . . ) is the vector of the remaining Qj-s, in non- 
increasing order. We denote by /i the distribution of C = (Co, Ci, C2, . . . ). 

Theorem 2. The distribution /i is invariant under the above split- and-merge trans- 
formation. 

Proof: By definition a random partition C with distribution fi can be considered 
as follows. Let Wi,W2, ■ ■ ■ be i.i.d. uniform random variables on [0,1] as in the 
definition of PD(1). Because the size biased permutation of PD(1) is GEM(l), we 
can suppose that for the active component Co = Wi holds and (Ci, C2, . . . ) is the 
decreasing rearrangement of {{1 — Wi)W2, {I — Wi){l — W2)W3, . . .) . Let be the 
distribution of the random partition obtained by the appHcation of a stirring step 
to C. 

If [/ < Wi for the [0, 1] -uniform random variable U, then the new non-active 
components are {Wi — U,{1 — Wi)W2, (1 — W^i)(l — VF2)W3, . . . ) in decreasing order. 
Conditionally on {U < Wi} and on U, the variable Wi is uniform on [U, 1], thus 
the vector of the non- active components has PD(1) distribution scaled by {1 ~ U). 
It yields that v conditioned on {U < Wi} and on U is the same as /i conditioned 
on the active component having size U. 

If U > Wi, then a coagulation occurs. Conditioned on {U > Wi} and on the 
value of Wi , the size of the component which merges Co has uniform distribution 
on [0, 1 — Wi], because it is a size biased component. We get the same distribution, 
if we choose this component merging Cq to be of length U — Wi. Conditionally on 
{U > Wi} and on U the rest has PD(1) distribution scaled by (1 — U). Thus, a 
sample from ly conditioned on {U > Wi} and on U has an active coordinate of size 
U and the remaining components with a scaled PD(1) distribution. 

Hence, a vector with distribution v can be obtained by sampHng U uniformly 
on [0, 1], taking the active coordinate of length U and taking a scaled PD(1) distri- 
bution on the rest. It shows that — fi, as required. 

Theorem 2 proves that /i is a stationary measure for our process, but it is not 
at all clear if this is the unique stationary measure. The proof of this would be the 
analogue of Schramm's result in pUj. 

A possible generaHzation of the model studied in this paper is the multiple stir- 
ring. It means that we consider more than one stirring particles. For a fixed number 
k of stirring elements an analogous limit theorem can be proved with a coupling 
similarly to Theorem 1. The case, if the number of the stirring elements depends 
on the size of the set [n], might also be worth studying (for example with k{n) = n" 
where < a < 1). Of course, we need different scaling of time and space in this 
case. 
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An open question is for our original model to establish after how much time a 
permutation can be regarded as a random permutation chosen with uniform dis- 
tribution, if it can be regarded at all. The solution of the problem in this simply 
describable model is not obvious in the least. For more about this problem in similar 
models see [3j. 
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